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Using a generalized least square procedure, a method is developed to estimate diffusion coeffi- 
cients, Dij, from concentration measurements made in a set of experiments of varying duration but with 
the same initial conditions. The method requires neither approximations nor restrictions in the flow 
equations. The Dij and their errors are calculated simultaneously with the adjustment of the weighted 
observables — concentrations of the two solute components and time. A procedure is described to 
insure convergence in a given experimental set. The effect of the duration of the longest experiment, 
and errors in the observables on the estimated values of D tj and their errors, are investigated. Applica- 
tion of the method to available data yielded values for D tj in agreement with those previously calculated 
by more involved procedures. 
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1. Introduction 

The methods reported for the calculation of diffusion 
coefficients in ternary systems from diaphragm cell 
data require either restrictions in the experimental 
conditions [1, 2] 1 or approximations in the flux 
equations [3]. A new method is reported here that 
obviates the limitations of previous procedures when 
concentrations of the diffusing solutes are directly 
measured in the cell compartments. Although the 
present method was derived independently of a pro- 
cedure recently reported by Cussler and Dunlop [4], 
the two approaches to the problem are similar. The 
present derivations, however, allow a more complete 
and direct treatment of the subject, particularly for 
the calculation of errors. 

2. Theory 

The general equations that define diffusion flows, 
(J/)i, in ternary systems relative to a volume-fixed 
frame of reference are [5, 6]. 
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2 fif. 

(Ji)v — Y,(P»k-z? (*=1,2) 



(1) 



in which (Z)y)v are practical diffusion coefficients and 
Cj is the concentration of component j. Integration of 
eq (1), considering the conditions for the diaphragm 
cell, assuming linear concentration gradients and no 
changes in volume during the diffusion process, gives 
the general solution [7] 



ACi = A i e-v»h< + B l e-P m 2< (i=l, 2) 



(2) 



in which AC; is the concentration difference (molar 
scale is used throughout these derivations) of com- 
ponent i between the two compartments at any time t, 
f3 is the cell constant and the other quantities are 
defined by 



A t = 



B,= 



AC°(D ii -m 2 )+D ij AC° 
rri\ — rriz 

Aq(D ii - mi )+D i Aq 

m-2 — m\ 



i=l, 2 (3) 
7=1,2 

i+j (4) 



s[Z>n + £>22 + V(D„ + Z) 22 ) 2 -4(Dn£> 22 -D 12 D 21 )] 



m 2 = |[^ 11 + £' 22 -V(D ll + D 22 ) 2 -4(Z) 11 £>22-£> 1 2Z)2i)] 

(5) 



43 



In eqs (3) and (4), AC° denotes concentration differ- 
ences at t = 0, and Dij are the integral cell diffusion 
coefficients. 2 

The calculation of Dij is accomplished by a general- 
ized least squares procedure [8] using eq (2) as the 
mathematical model. Thus, the condition equations 
are of the form 

F= (Zi-Z 2 ) -A l e-P m i T -B l e-f> m 2 T (6) 

L= (7i-y 2 ) -Ate-fif-Btfi-**** (7) 

in which X and Y represent the concentrations of 
components 1 and 2, respectively, in the compart- 
ments indicated by their subscripts at the observation 
time T. 3 

If observations (determinations of X and Y) are 
made at the end of N diffusion experiments (with the 
same initial conditions) of duration 7\, T 2 , . . ., 
Tn-, then there will be 2N condition equations. The 
observations corresponding to different values of T 
are statistically independent, since for each duration, 
a new experiment is conducted. The number of condi- 
tion equations may not be a multiple of 2, however, 
if either X or Y is not determined at the end of some of 
the experiments. 

When observations are made in only one of the 
compartments, for example compartment 2, the 
condition equations take the form 

F'=(AX +2Xl-2X 2 )-Aie-e m i T -B i e-f im 2 T (8) 

L' = (AY° + 2Y° 2 -2Y 2 ) -A 2 e-e m i T -B 2 e-V m 2T (9) 

in which AX° and AY° are concentration differences for 
the two components at time T= 0. 

The aim of the procedure is to compute the adjust- 
ments of the observables and the corrections for the 
initial estimates of the diffusion coefficients in such a 
way that 

F^(x[ n \ x 2 n \ t< n \ Ax°, Ay°, D lu D 12 , D 2U D 22 ) -0 

(n = l,2, . . ., N) (10) 
and 

/>>(;#>, y£\ |« A* , Ay°, D lu Z) 12 , D 2U D 22 ) -0 (11) 

in which x, y, t, A#°, and Ay° are adjusted observables, 
and Dij are the final estimates of coefficients. This 
adjustment is done subject to the least squares cri- 
terion that the weighted sum of squares of the residuals 
(observed value — adjusted value) of the observables 



2 Usually, diffusion experiments in a diaphragm cell will yield values for integral diffusion 
coefficients. However, if AC° is selected sufficiently small, the calculated Du will approach 
the values for differential diffusion coefficients. The present treatment will apply equally 
well for such conditions. The functional relationships between the integral Dij and con- 
centration of the components in ternary systems are not known. Attempts have been made 
to solve this problem for some particular systems [P. F. Mijnlieff and H. A. Vreedenberg, 
J. Fhys. Chem. 70, 2158 (1966)] but a general analytical solution involves formidable 
mathematical difficulties. 

3 After the preparation of the present manuscript, J. Mandel, National Bureau of Stand- 
ards suggested the use of the concentration differences in eqs (6) and (7) as observables, 
instead of working with individual concentration terms. The present procedure would apply 
equally well for this case. 



X, Y, T, AZ°, and AF° be a minimum. In these deriva- 
tions, the cell constant, /3, is considered as an "external 
parameter," that is, a quantity obtained from independ- 
ent experiments, and therefore, not susceptible of 
adjustment. However, the error in /3 is taken into 
account in the estimation of the error bounds for 
Dij, as explained below. 

Linearization of the condition equations (6) and (7) 
is done about the observed points (Xi,X 2 , 7\ AZ°, AF°, 
D° u , D° 12 , D° 2U D° 22 ) and (F„ Y 2 , T, AX°, AF°, DJ U D° l2 , 
D 2X , D22) using the best available initial estimates, 
D°j, for the diffusion coefficients. A procedure to obtain 
good initial estimates for Dij is described later. 

The 2N reduced condition equations, resulting from 
the linearization mentioned above, may be written, 
in matrix notation, 



BV+Pa = M 



(12) 



in which B is the matrix (2W X [57V + 2]) of the deriva- 
tives of the condition functions [6] and [7] with 
respect to the observables, V is the column vector 
( [57V + 2] X 1) of the residuals, P is the matrix (27V X 4) 
of the derivatives of the condition functions with 
respect to Dij, a is the column vector (4 X 1) of the 
corrections for the Z)°- (i.e., a tj = D°. — Dij) and M is the 
column vector (2/VX 1) of the values of the condition 
functions. The condition functions and their deriva- 
tives are evaluated using the experimental values of 
the observables and the initial estimates for the dif- 
fusion coefficients. 

The sum of squares of the residuals, S 2 , is given by 



S 2 =VWV 



(13) 



in which W is the diagonal matrix of the weights, W^s, 
for the 57V + 2 observables and the tilde represents 
transposition. These weights are calculated from 



w. 



obs " 



_<L 



a; 



(14) 



obs 



in which cr obs is the standard deviation for the ob- 
servable in question and cr int is an arbitrary constant 
conveniently selected for calculation purposes. 

Minimization of S 2 (eq (13)) subject to the conditions 
(12), effected by the use of a Lagrange multiplier, 
A, and the subsequent elimination of V, results in 
27V + 4 normal equations with 2N unknown X's (the 
elements of column vector (A), and four unknown a { j. 
As shown by Guest [9] this leads to the following 
equations: 



[J 



R P 

P 



A 

Q 



M 





or 



(15) 
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in which R = BW^B, of dimension 2N X 2iV. 
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The solution of the above system of equations yields 



(16) 



Since G~ J is symmetric, like G, it can be partitioned 
similarly as follows: 



A _ 


= G-> 


M - 


a 








4 2 



Q H 

C 



(17) 



Then, the solutions for A and a can be written as, 

A = QM 

a = HM 

The values thus obtained for atjj can be now used to 
calculate new estimates for Dij, since 



D«=D?,-a«. 



(18) 



These new estimates of Dij are used iteratively as 
initial estimates (substituting for Z)°) until convergence 
in S 2 , calculated by eq (13), is obtained. At this point 
the final residuals are calculated from V=W~ X BK. 
The error vector A, (ADn, AZ), 2 , AZ) 2 i, AZ) 22 ), can be 
found from the 4X4 error submatrix, E=—C, of G~ l : 



A = < 



(19) 



in which the standard error, cr ext , is given by 
V{S 2 /(2yV — 4)}, 2/V-4 being the degrees of freedom. 

A discrepancy between the ratio cr 2 xt /o-f nt and unity 
is an indication of disagreement between the variances 
assumed in the calculations (based on previous ex- 
perience) and the error variances actually existing in 
the experiment [10]. 

The errors as obtained from eq (19) do not reflect 
the influence of the error, cr^, in the independently 
determined cell constant /3. Therefore, calculations 
are repeated using successively the values /3 ■+■ crp 
and /3 — ap for the cell constant. A t-test is then per- 
formed between the Dij obtained in this fashion and 
those obtained when (3 was considered free of error. 
If the Mest indicates no significant difference, at a 
selected level of significance, the errors obtained 
for the Dij with the assumption that crp = may be 
adopted. Otherwise, a confidence band for the Dij 
should be estimated. 

A requirement for the convergence of the iterative 
procedure in the calculation of Dij is that the initial 
estimates of these parameters be close to the final 
estimates. Various procedures have been suggested 
[11, 12, 13J to insure convergence in problems of this 
nature. In this investigation it was found that, for 
ternary systems, good initial estimates for only two 
of the four coefficients are required to insure con- 
vergence, their identity being determined by the initial 
experimental conditions. Three cases of initial con- 
ditions are possible in a ternary system: (1) Both dif- 



fusing components have nonzero concentration gradi- 
ents; (2) Component 1 has a nonzero gradient and 
component 2 has a zero gradient; (3) Component 2 has 
a nonzero gradient and component 1 has a zero 
gradient. The two diffusion coefficients for which good 
initial estimates are required in each one of the three 
cases are: for case (1) D u and Z) 22 ; for case (2) D n 
and D 2 u and for case (3) Z) 22 and Z)i 2 . The initial esti- 
mates for the remaining 2 coefficients in each case 
can be taken as zero. In the flux equations (1) the main 
contribution to the flux of one component is given by 
those terms involving the component for which there 
is an initial concentration gradient. Then, ignoring, as 
a first approximation, those terms that contribute the 
least to the diffusion fluxes, performing the corre- 
sponding integrations, the resulting expressions may 
be written in the logarithmic form, 



i AC, 
ln-££5=-DvPt 



i=l, 2 

;=l,2 



(20) 



which covers the three cases given above. The required 
initial estimates of the diffusion coefficients are then 
obtained by a simple linear regression using the ap- 
propriate experimental data. In this way, convergence 
was obtained after 4 to 6 cycles in all the cases tested. 



3. Test of the Method 

The procedure was thoroughly tested by construct- 
ing artificial "experimental" data. The condition 
functions (8) and (9) were used assuming a value of 
zero for X% and Y%. Exact solutions of eqs (2) were 
obtained for different arbitrary sets of D#, /3, and 
initial concentration differences. The "true" points 
thus obtained were transformed into "experimental" 
observations by the use of a table of normal deviates, 
and the foregoing method was applied to these data. 4 
Thus the following factors affecting the estimated 
values for Dij and their errors were investigated: 
(a) duration of the longest experiment in a given set 
of experiments, (b) number of experiments (total num- 
ber of observations) in sets having the same duration 
for the longest experiment, and (c) errors in the 
observables. 

For the case of binary systems, expressions have 
been reported [14, 15] to calculate an optimum ex- 
perimental time at which the fractional error in the 
diffusion coefficient becomes minimum. Derivation of 
similar expressions for ternary systems involves such 
complexity that it is practically unmanageable. In 
the present investigation this subject was studied by 
calculations of Dij for different sets of experiments, 
varying the duration of the longest experiment in 
each set. A sample of the results obtained is shown 



4 All the values given in the tables in this paper were calculated assuming that the rela- 
tive precision in the determination of the observables (except time) remained constant. 
This situation is often found when concentrations are directly determined by chemical 
analyses; the volume of the aliquots is selected so as to analyze essentially the same amount 
of component in the samples. The adjustment procedure, however, is general and satis- 
factory results are also obtained by considering absolute precision. 
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in table 1. The time given at the top of each column 
corresponds to the longest experiment in a set con- 
sisting of eight experiments, the time increment for 
each successive experiment being one-eighth of the 
longest experiment. Values for the percent decrease 
of the initial concentration differences of components 
1 and 2, Ri and R2, at the end of the longest experiment 
for a particular set are given in the second and third 
rows. In order to facilitate the interpretation of these 
calculations, the same random numbers were used to 
estimate the "experimental" observations for each set. 
For this reason, the ratio cr| xt /crf nt is fairly constant 
in this table and, as might have been anticipated, its 
value is close to unity. 

Table 1 shows that as the duration of the longest 
experiment increases, the estimated Dij approach the 
values for the "true" diffusion coefficients. Further- 
more, the last row in the table shows that there is a 
minimum in the sum of squares of the errors in Dy, 
SAfj, corresponding to a set with an experiment of 
maximum duration between 200 and 600 hr. This type 
of calculation was also performed with a number of 
systems varying the initial conditions and the (true) 
values for D X j. In all cases a minimum was found but 
its position varied from system to system and there 
was no apparent trend correlating the minimum with 
the parameters used or calculated. In some of the 
systems studied, the approach to and departure from 
the minimum value of 2A| occurred very gradually 
as the duration of the longest experiment increased. 
In other systems, however, the valley region was very 
narrow. Typical examples of these features are shown 
in figure 1 in which the quantity £A| is plotted against 
the duration of the longest experiment in various ex- 
perimental sets. Initial conditions and related param- 
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Figure 1. Effect of duration of longest experiment on errors in Dy. 

For both plots: Number of experiments per set, initial conditions, errors in observables, 
and cell constant, as given in table 1. For la, "true" Dij as in table 1. For lb, "true" 
D u x 10 s : D n = S.0, D 12 = 1.0, D 2l =-2.0, D 22 = 4.0. 



Table 1. 


Effect of duration of the longest experiment on the 


estimated values and errors of the diffusion coefficients 


Duration of 


24 


120 


216 


312 


408 


504 


600 


800 


True diff. 


longest exp., hr* 


















coeffs. 


/^, = 100 AX/AX° 


83.2 


39.5 


18.3 


8.1 


3.4 


1.2 


0.2 


-0.3 




R 2 =100 AF/AF° 


90.5 


61.2 


42.0 


29.1 


20.4 


14.3 


10.1 


4.9 




D n x 10 5 


1.576 
±0.155 


1.935 
±0.033 


1.959 
±0.021 


1.967 
±0.018 


1.971 
±0.016 


1.973 
±0.016 


1.974 
±0.017 


1.975 
±0.019 


2.000 


D 12 x 10 5 


0.888 


0.296 


0.255 


0.240 


0.233 


0.229 


0.226 


0.221 




±0.252 


±0.049 


±0.030 


±0.023 


±0.021 


±0.020 


±0.020 


±0.022 


0.200 


D, x x 10 5 


0.062 


0.103 


0.103 


0.103 


0.103 


0.103 


0.103 


0.103 




±0.052 


±0.011 


±0.008 


±0.006 


±0.006 


±0.006 


±0.006 


±0.006 


0.100 


D 22 x 10 5 


1.062 
±0.085 


0.997 
±0.017 


0.997 
±0.011 


0.998 
±0.008 


0.998 
±0.007 


0.998 
± 0.007 


0.998 
±0.007 


0.999 
±0.007 


1.000 


O-ext 


0.72 
97.614 


1.07 
3.909 


1.08 
1.489 


1.08 
0.946 


1.08 
0.770 


1.07 
0.724 


1.08 
0.744 


1.09 
0.916 




°fnt 

(IAD x 10« 







Initial conditions: AX° = 0.025; AF°=0.015. Cell constant: = 0.1. 
Total number of experiments in a set: 8. 

Adopted errors: £ A , >=±0.2%; £ r = ±60 s; E&*, ay° = ±0.1%; #0=0.1%. 
*Same set of random numbers used for each column. 
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Table 2. Effect of number of experiments in a set on the estimated values and errors of the diffusion coefficients* 



Number of 

experiments 

in a set 


3 


5 


8 


True diff. 
coeffs. 


D n XlO s 
D12XJO 5 

D, x X 10 5 

D 12 x 10 5 


2.005 
±0.027 

0.197 
±0.029 

0.114 
±0.009 

0.986 
±0.010 

0.58 


1.938 
±0.026 

0.261 
±0.028 

0.088 
±0.009 

1.012 
±0.010 

0.58 


2.051 
±0.015 

0.151 
±0.018 

0.097 
±0.005 

1.003 
±0.006 

0.47 


2.010 
±0.022 

0.202 
±0.025 

0.096 
±0.008 

1.004 
±0.009 

0.99 


2.002 
±0.010 

0.189 
±0.012 

0.103 
±0.003 

0.994 
±0.004 

0.40 


1.973 
±0.016 

0.228 
±0.020 

0.103 
±0.006 

0.998 
±0.007 

1.08 


2.000 
0.200 
0.100 
1.000 


<rL 





Initial conditions: A*° = 0.025; AF° = 0.015. Cell constant; 0=0.1. 

Longest experiment in each set: 516 hr. 

Adopted errors: £ A , y = ±0.2%; £Y = ±60 s; E hx „. M - ( , = ±0.1%: E = ±0.1%. 

^Entries in subcolumns are extreme values of Dn obtained with 6 to 10 sets of randomized data. 



eters are given in the figure caption. The lower curve 
corresponds to the example given in table 1; however, 
both plots are restricted to the regions adjacent to the 
minimum. Values for Dij estimated from the experi- 
mental sets represented by the points between arrows 
in figure 1 do not differ appreciably. 

The actual position of the minimum in £A^. for a 
given system shifts for different values of the cell 
constant, /3, but it always occurs at approximately 
the same values of R\ and R 2 . In the absence of 
analytical expressions to calculate optimum diffusion 
times in ternary systems, the procedure reported here 
can be used to determine the duration for which a 
minimum in the errors for Dij is obtained. This ap- 
proach, of course, requires estimates for the values 
of Djj, errors in the observables, initial conditions, 
and the cell constant. With these parameters, curves 
like those in figure 1 can be computed and a decision 
about the diffusion time can be made. 

The effect of the number of experiments in a set 
on the estimated values of Dij is shown in table 2. 
Calculations for this table were made using the same 
initial conditions and set of parameters used in table 
1. The longest time used in any particular set was 
516 hr, and the time increment for the TV singular 
experiments was 516/N. The calculations were carried 
out several times for each set using different sets of 
random numbers for the estimation of the "experi- 
mental" observations. Results for the two sets of 
calculations displaying the maximum differences in 
Dij are given in the table. It is apparent that as the 
number of experiments in each set is increased, the 
variations in Dij and A/j due to experimental errors 
decrease although the duration of the longest experi- 
ment in all the sets was the same, 516 hr. 

The errors in the observables not only affect the 
estimated values for Dij and A,j, but also the time at 
which a minimum in XA'f, is obtained for a given 



system. Such effects are shown in table 3; the values 
for Dij were calculated from sets having an experi- 
ment with a maximum duration equal to the time for 
which a minimum in SAfj was observed. With large 
errors in the observables, the values for Dij obtained 
by the adjustment procedure may be very inaccurate; 
this situation is particularly true for the cross-term 
diffusion coefficients. Although definite shifts in the 
location of the minima in SAfj are apparent, they are 
hardly significant for experimental planning. 



Table 3. Effect of errors in observables on the estimated values 
and errors of the diffusion coefficients 



E X , y 


2.0%* 


0.2%* 


0.02%* 


True diff. 
coeffs. 


Z) U X10 5 
D 12 X10 5 

D21 x 10 5 

D,> x 10 r > 
(SAfJxlO 13 


1.742 
±0.147 

0.465 
±0.181 

0.130 
±0.055 

0.978 
±0.068 

62.025 
532 


1.973 
±0.016 

0.228 
±0.020 

0.103 
±0.006 

0.998 
±0.007 

0.723 
516 


1.997 
±0.002 

0.203 
±0.002 

0.100 
±0.001 

1.000 
±0.001 

0.008 
536 


2.000 
0.200 

0.100 
1.000 


Duration of 
longest experi- 
ment, hr 







Initial conditions: AZ° = 0.025; AF° = 0.015. Cell constant: = 0.1. 
Total number of experiments in a set: 8. 
Adopted errors: £ax°,ay° = ±0.1%; £^ = ±0.1%; £ 7 =±60 s. 
*Same set of random numbers used for each column. 
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The adjustment procedure described in this paper 
was applied to the diffusion data reported by Kelly 
[1] for the ternary system mannitol(l) — NaCl (2) — H 2 0, 
using the condition eqs (6) and (7). Kelly had ob- 
tained the values for Dij from the results of two 
experimental arrangements with initial conditions 
given at the head of table 4. His results and the 
results obtained by the present procedure are shown, 
side by side, in the same table. 

There is satisfactory agreement between the 
results obtained by the present procedure and those 
obtained by Kelly for systems where he could cal- 
ulate diffusion coefficients, with the exception of 
the value for D u in the system 2— A. The values for 
Dij estimated from Kelly's experimental data for 
system 1— A should have agreed, in theory, with those 
obtained from the experimental data for the system 
1-B. The cause of these discrepancies is that, usually, 
when one of the components has zero initial concentra- 
tion gradient, the amount diffused is so small that a 
reliable estimation of the Dy is impossible. 

The discrepancies in systems 2-A and 2-B are 
explained on the same basis; however, the large 
value for the ratio cr| xt /cr? nt for system 2-A suggests 
that the error variances used for the observables are 
very different from those actually existing in this 
system; for this reason, the values for Dij obtained 
from this experiment are questionable. 

Cussler and Dunlop [4] also estimated Dij from 
Kelly's data. It is not clear in their publication whether 



they combined the data for systems 1-A and 1-B or 
estimated Dij from each arrangement separately. The 
values they obtained are in good agreement with those 
obtained by the present method and reported in table 
4, with the exception of the value for D 2 i which is 
about 5 percent lower than the figure reported here. 
Significant differences, however, do exist between 
the errors reported by Cussler and Dunlop and those 
obtained in the present investigation. This is not 
surprising because in their unweighted least squares 
procedure, those investigators arrived at estimates of 
errors for the Dy by a bound method. Their reported 
error figures correspond to the maximum difference 
between the Dij calculated considering the nominal 
value of the observables and the Dij obtained through 
similar calculation but using the values of the ob- 
servables as the nominal values minus their errors 
and plus their errors, respectively. The procedure 
given in this paper to calculate A,j is not only time 
saving but also it has a sound statistical basis. 



The authors are indebted to L. J. Gosting, Institute 
for Enzyme Research, Wisconsin University, Madison, 
Wisconsin, for making available to them the material 
of reference (1) and to R. P. Wendt, Department of 
Chemistry, Loyola University, New Orleans, Louisiana, 
for his suggestions and critical review of the 
manuscript. 



Table 4. Comparison of Kelly's estimated values of the diffusion coefficients with those estimated by the present procedure using Kelly's 

experimental data for the system mannitol (X) — NaCl (Y) — H 2 



Average 
cone. 


X=0.2M; F=1.5A/ 


X=*0.2M 


; Y=3.0M 




Initial 

conditions 


&X° = 
AF° = 


-0.4 
= 0.0 


AF° = 


= 0.0 

= 0.5 


A* = 0.4 
AF° = 0.0 


&x°-- 

AF°: 


= 0.0 
= 1.0 


System 


1 -A 
5 experiments 


1 -B 

4 experiments 


2-A 
6 experiments 


2-B 
3 experiments 




This 
paper 


Ref. 

(1) 


This 
paper 


Ref. 

(1) 


This Ref. 

paper (1) 


This 
papei 


Ref. 

(1) 


D„xl0 5 


0.557 
±0.003 


0.559 
±0.004 


1.148 
±0.246 




1.006 0.448 
±0.016 ±0.006 


0.439 
±0.625 •' 










D V2 x 10 5 


0.025 
±0.039 




0.018 
±0.002 


0.016 
±0.010 


0.512 


0.013 
±0.003 


0.013 






±0.066 


± 0.008 


Z)„ X 10 5 


0.232 
±0.017 


0.232 
±0.001 


0.138 
±1.918 




0.494 0.456 
±0.072 ±0.003 


-0.040 
±9.969 










D„xl0 5 


1.556 
±0.230 

5.20 




1.388 
±0.013 

4.32 


1.390 
±0.002 


1.122 
±0.291 


1.469 
±0.048 


1.469 


°-'L 




± 0.003 








167.18 


38.98 




<t 











Ejr,r> = ±0.02%;E x , 



:0.02%;£Y=±60s;£tf = 0.1%. 

48 



4. References pp - 128 ~ 147 (J° hn Wile y & Sons ' Inc - New York ' New York - 

1943 ) (Reprinted by Dover Publications, New York, New 

[1] Kelly, F. J., Ph. D. Thesis, University of New England, Armidale, rn , r . i* r> kt i \/i *u j en l-»*- «j-7n qoi 

NSW A ftl'i (19611 ' ' Guest, P. G., Numerical Methods of Curve rating, pp. 370-381 

roi v ii" "f M i j c. i u u li i * ru: inn u (Cambridge, University Press, Cambridge, England, 1961). 

2 Kelly, r. J., and Stokes, R. H., Electrolytes, pp. 96-100, B. Mn , D . D ^ du D • .n on-7 /moo. 

Pesce Ed (Peraamon Press I td Oxford 1962 [10J BlFge ' R ' T " PhyS> Review 40 ' 207 < 1932 )' 

rol „ , \ ¥ V JT u i J i 'ni / u A oni:/in^ [HJ Box, G. E. P., Bull. Inst. Internat. Statistique 36,215(1958). 

[3 Burchard, J. K., and Toor, H. L., J. Phys. (.hem. 66,2015 1962). \ 10 \ u , u r* t u o o*n/irw:i\ 

L n , r t i j rk i d i i du Vl t^ 1121 Hartley, H. O., lechnometncs 3,269(1961). 

14 (.ussier, E. L., Jr., and Dunlop, P. J., J. Phys. Chem. 70, moi c* j t r- V ui n a j d u da i r-u di 

iQQnnoAA\ [13] Strand, T. G., Kohl, D. A., and Bonham, R. A., J. Chem. Phys. 

188UUVOO). oq iQ07nQ6^^ 

[51 Onsager L Ann NY Acad Sci 46 241(1945) ' u^uoj. 

L i^-i j t r- d u • d t 'rk 'i d t r- . t t [14-1 Robinson, R. L., Jr., Edmister, W. C, and Dullien, F. A. L., 

[6| Kirkwood, J. G., Baldwin, K. L., Uunlop, K. J., Costing, L. J., j p, p. 69 258H965^ 

m rin d Wp G R J ^ hem R ?h r 33 J r ! 5 /, 96 ^ n I r I t 15 l Va « Git, A. L," and Adamson, A. W., J. Phys. Chem. 68, 238 

[7 1 Gilhland, E. R., Baddour, R. F., and Goldstein, D. J., Can. J. L J /-iQfid* 

Chem. Eng. 35, 10 (1957). < 1W * p _ Q . 

[8] Deming, W. E., Statistical Adjustment of Data, pp. 49-58 and (raper MAI 566) 



49 

324 194 0-69— 4 



